version 30.9.2012 



Solution of the relativistic bound state problem for hadrons 
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Abstract 

To overcome the infrared divergency problem of the Standard Model of parti- 
cle physics, a second order Lagrangian - containing two boson fields - has been 
evaluated, which includes all mechanisms necessary to form a relativistic bound 
state. Strikingly, in the application to hadrons an overall consistent description 
of qq flavour systems is obtained only by assuming massless elementary fermions 
(quantons). Further, the important but not well understood confinement potential 
is obtained in an Abelian theory, indicating that confinement is not due to colour 
but is a general property of relativistic bound states. 

Two mesonic systems are discussed, a;(782) and charmonium J/*I'(3097) includ- 
ing excited states, for which a satisfactory account of the absolute masses and a 
qualitative understanding of their decay widths is obtained. 

PACS/ keywords: 11.15.-q, 12.40.-y, 14.40. Cs/ Bound state description of hadrons 
based on a second order Lagrangian with massless fermions (quantons) and two- 
boson coupling. Confinement and boson-exchange potential, masses and widths of 
a;(782) and J/^'(3097). No coupling to scalar Higgs fields needed. 

To study fundamental forces, nature provides us with hadrons and leptons, which form the 
constituents of matter, but also with composite systems, nuclei, atoms, and gravitational 
states in form of solar and galactic systems. For the description of these stable and massive 
systems as bound states a relativistic theory is needed, since the elementary constituents 
of these states are relativistic. However, relativistic bound state problems are generally 
difficult to solve, see e.g. Salpeter and Bethe [1], and could not be tackled so far for 
particle bound states (see the discussion in ref. [2]). 
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Instead, powerful effective theories have been developed, which are contained (except 
gravitation) in the Standard Model of particle physics [3J (SM). In this framework particle 
bound states are included effectively by a number of parameters including nine masses 
of 'elementary' particles. Many experimental data are described, but the underlying 
mechanisms leading to particle bound states are hidden. Attempts to solve the relativistic 
bound state problem for hadrons in the framework of the SM were not successful [2] due to 
the infrared divergency of quantum chromo dynamics (QCD), which makes it impossible 
to derive hadrons of finite radius. 

So far, hadronic bound states could be simulated only in the lattice approach ^ . However, 
in this complex empirical approach the mechanisms which lead to bound states cannot 
be clearly identified. As an example, the confinement potential - of prime importance for 
hadron structure - is produced in lattice calculations, see e.g. the review in ref . [5] , but the 
underlying mechanisms are still unclear, e.g. whether it can be understood in the picture 
of colour strings or some complicated effect of gluon binding. For a better understanding 
of fundamental forces it is therefore imperative to study all possible mechanisms, by which 
bound states of relativistic particles can be formed. This should lead to a finite theory, 
in which the infrared divergency problem of the SM is removed. 

Different from bound state potentials in a non-relativistic theory, those of relativistic 
particles can be created only in complex processes, involving an interplay of fermions and 
bosons. By studying possible mechanisms, which could lead to such bound states, it was 
found that this appears possible in qq creation during the overlap of two boson fields. But 
this is not sufficient, it has to make sure that this system does not separate and forms a 
two-boson/two- fermion bound state. 

To describe these processes, a second order SM Lagragian has been used, extended by two 
boson fields. Formally, such a Lagrangian yields solutions only, if the time component of 
both boson fields is the same. This allows to replace the relativistic fermion four-vectors 
by corresponding ones in three dimensional space, leading naturally to a non-relativistic 
system, for which bound state potentials can be deduced. The dynamics of such a system 
appears to be rather complex, but bound state properties can be studied more easily. 

To facilitate such a description, instead of an extension of the non-Abelian QCD La- 
grangian a simpler Lagrangian of Abelian form was used. This may to be appropriate. 
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since hadron systems are colour neutral. However, such a description should fail com- 
pletely, if the confinement of hadrons could be generated only by the colour structure of 
QCD. 

The Lagrangian has been used in the form 

C = ^^ ij^D^^D^D"^ - V^.F'^^ , (1) 

where m is a mass parameter and \E' a two-component fermion field \E' = (\E'"*' \E'°) and 
^ = (\E'~ with charged and neutral part. Vector boson fields are contained in 
the covariant derivatives = — igA^ and the Abelian field strength tensor F^^ = 
Qfj-j^f _ Only couplings to the charge {g = gc) and spin {g = gs) of fermions are 

considerecl^. In the bound state problem we deal with a finite quantum field theory, so 
renormalisation is not needed. 

We insert Df" = d^" - igA^ and D^D"" = d^d"" - ig^A^d" + d^A'') - g'^A^A^ in eq. (HD and 
obtain for the first term of C 

^z 

■ 2 

m"^ 

• 2 '2 3 

- j^A'^d^A-'^ - j^d^'A.A-'^ - 7^AM,A^^ . (2) 

The gauge condition d^A^ = used for simpler Lagrangians (as in QED) is replaced in 
our case by dld^A") = 0. Also for the fermion field we have to require a cancellation 
of the second and higher derivatives. By this condition unphysical solutions are avoided, 
which can arise [6] in Lagrangians, which contain higher derivatives of the fermion field. 

In eq. (|2]) the number of field derivatives and boson couplings varies between the first and 
last term. This shows that the various terms are related to different kinetic situations, 
pointing to a rather complex dynamics of the system. 

Contributions to stationary solutions have been studied by evaluating fermion matrix 
elements (or ground state expectation values) of the field operators 0^(g) in eq. (|2]). 
These can be written by J^{p' —p) = < g-s.\ K{q) \g.s. >~ i^ip') K{q) 4'{p), where ipip) 
is a fermionic wave function ipip) = ^^^^(^1)^(^2) and K{q = p' — p) a kernel, which is 
^in the application discussed below only charge coupling is needed. 
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expressed by K{q) = ^2{n+i) [0]^{(l)gfipOp{q)], see e.g. [7], where n is the number of boson 
fields and derivatives in eq. ([2]) (in the present case n=3) and (7^p the metric tensor. 

In the following the essential points of the derivation of matrix elements and bound state 
potentials are given. For the construction of stationary states we expect contributions 
mainly from terms of the Lagrangian ([2]), which contain static fields (without derivatives). 
This is the case only for the last term Ci^s = g^'JfiA^A^A'^'^ and leads to a matrix 

element Aisg, which contains three boson fields on the right and left 

-a^ - 

Msg = ^{v')l,A>^{q)AM)A''{q) g,, AMA''{q)l,A''{q)^{v) , (3) 

where a = g'^/4:7c. A comparable matrix element in the (first order) SM may be written 
in the form AisM = ^ V'(p')7m^^(Q') S'mp 1pA^{q)'ip{p), giving rise to a (boson-exchange) 
interaction of vector structure v^{q) ~ aA^{q)AP{q) (but only for equal times of the 
two boson fields, which means in the non-relativistic limit). Differently, in eq. (|3]) three 
interactions are involved, one of vector structure Vv{q) ~ aA^{q)A'^{q) or ~ aA,y{q)AP{q), 
the other two with a coupling of two boson fields to a scalar (~ aA,j{q)A'^{q) and ~ 
aA^{q)A''{q)) or to a vector (~ «A^(g)A^(g) and ~ aA^{q)Af'{q)). 

In a dual picture the two-boson fields, which appear twice (on the left and right side of 
Aisg) can be regarded as interactions, but also (analoguous to the fermion wave function 
ip{p)) asbosonic (quasi) wave f unction Jf|ws(g) = 4-ylj^(g)A''(g) andwv{q) = ^Ap(q)A''(q). 
However, this formal construction does not lead to a boson density for excited systems, 
which can be interpreted as probability. The physical picture of M.3g is that for the 
fundamental state of a relativistic bound state system the fermions interact only inside the 
two-boson density (which leads to boundary conditions discussed below) and feel therefore 
three interactions. A single boson-exchange interaction is possible only in dynamical 
situations, see the terms 2-4 in eq. ([2]), which do not lead to a bound state potential. 

The 7-matrices as well be removed (contracted) by adding a matrix element 

with interchanged yU and p (according to |(7^7p + Iplp) = g^p)- Further, an equal time 
requirement of the two-boson fields (to reach overlap) allows to replace all fermion four- 
vectorcl by three-vectors in momentum or r-space. Correspondingly, the boson vectors 



^leading to boson (quasi) densities u>^(g) with dimension [GeT^]^ 
*in a {t, r) representation 
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are reduced to two dimensions. This yields 

— ct"^ - 

Msg = i^ip')Ws,vil) Ml) U!s,vil)^iP) ■ (4) 

The wave functions Ws{q) and w^{q) are given by the overlap of two boson fields with 
spin=0 and 1, respectively. The interaction Vy{q) has a structure similar to Wy{q). Writing 
the matrix element by A^3g = tplp') V^glq) ip{p), we obtain a three-boson potential 

VsTi^) = <M)vM ■ (5) 
Fourier transformation to r-space leads to a folding potential 

Vig\r) = / dr'wiy ) v^ir - r') . (6) 

Eq. ([6]) shows the typical form of a folding potential. Such a form has been used to 
describe elastic and inelastic hadron processes [8]. 

Concerning the fermion wave function %l){r) in the matrix element, here only il){r) ~ Ws{r) 
is assumed. In combined with the two potentials ([6]), this gives rise to two vector states, 
a 1~ fundamental state and one of significantly higher mass, which can decay into two 
mesons or baryons. Other wave functions are possible, e.g. a p-wave ^/'L=i(r) leading 
to 0"*" states, but these are not considered in the present paper. Since the fermionic 
and bosonic wave functions are finite, also the corresponding interactions are finite, in 
particular ^'^,(r). Consequently, solutions of the relativistic bound state problem exist in 
full space, as expected for a finite theory. 

To evaluate the potentials V"3*g^(r), the boson wave functions Ws^„(r) have to be determined. 
To achieve this, boundary conditions can be formulated by requiring that the interaction 
takes place inside the bound state volume of the g.s. As a consequence, the deduced 
potentials (E]) should be proportional to the density ip'^ij-) ~ wlir). For the potential 
V^g{r) this is not possible, but for the potential V^g{r) we can require 

ci wl{r) \V^^{r)\ , (7) 

if the density w1{r) has a radial node. Consequently, wlij-) must have a p-wave structure, 
which is translationally invariant and fulfills the condition < r^2 > =0. To achieve this, 
w^(r) is used in the form 

wlir) = wlir) $) = <(1 + /3R d/dr)wHr) ¥^,^0, $) , (8) 
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where wf,^ is determined from the normahsation / rdr w'^{r) = and (3R from the 
condition < r^„2 > = 0. 

Using the above constraints ([7]) and ([8]) this leads to a self-consistent boson wave function 
of the form 

Ws{r) = Ws, exp{-{r/hY} , (9) 

where Wg^ is determined from the normalisation An J rdr Wg{r) = fi^ and k from condi- 
tion ([7]). The extra interaction Vv{r) ~ Ws{r) is derived also from an expression similar to 
eq. (IE]), but here a translational invariance condition < r^, > = does not exist. Instead, 
(3R is adjusted to fit the boundary condition ([71) at large radii. 

To generate a stable bound state, the potentials V^g^{r) are not sufficient to keep the 
two boson fields confined. Inspecting the other terms of the Lagrangian ([2]), these show 
kinematic situations, in which bosons and/or fermions are in motion. Nevertheless, term 
6, which may be written in the form £i,6 = -'^^ -f^A'^^d^A")'^ - l^A^'Ay d""^, 
gives rise to another bound state potential. 

The first term of Ci^ leads a matrix element M.2gi which contains two boson fields and a 
derivative on the right and left 

With the boson gauge condition discussed above, we can write {dyA^{q)){daA"{q)) = 
^d^[da{A^A^Y]'^ . After contracting the 7-matrices and reducing the fermion and boson 
vectors by one dimension as discussed for A^sg, this yields 

M2, = ^^'(pO <iQ)d'w'M ^'{p) , (11) 

where ip'{q) and w'g{q) represent fermion and boson wave functions with a q-dependence 
different from ip{q) and Ws{q) in Ai^g. Since the two bosons are bound (momentarily), 
see eq. @, they have a kinetic energy d'^w'g{q)/2m. According to the virial theorem this 
implies also the existence of a static two-boson potential V2g{q). 

In a transformation to r-space the bosonic part of eq. ffTTj) gives rise to a Hamiltonian of 
the form 

a'^m<rl,2> ^d^w'^ir) 2dw'Xr)\ ,r f ^ n ^ r. > f ^ 
+ ~r~t^^ + '^^""^ '^''^ ^ ' '^''^ ' ^ ^ 
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where w'g{r) is the Fourier transform of w'^{q) given by Wg{r) = J qdq joiqr) w'g{q) and 
< > is the radius square of the boson density. The potential V2g{r) is given by 

where Eo is the energy of the lowest eigenstate. 

A similar potential involving w'^{q) deduced from Cij yields negligible contribution to the 
binding energy due to strong cancellations. All other terms of the Lagrangian ([2]) do not 
contribute to bound state potentials. 

Including the potential V2g{r) a final boundary condition for the matching of density and 
potentials can be formulated, requiring that the sum of all three potentials V^g^{r) and 
V2g(r) match the densities ip{r) ~ w^{r). However, V2g{r) cannot be added directly, since 
it is derived from a kinetic energy term, which itself is generated from the two-boson field 
Wg{r). Therefore, an adequate boundary condition can be formulated by 

C2 w',{r) = \Vigir)\ + |V3;(r)| + V^^ir) , (14) 

where V2g{r) ~ {a'^ /2m)w'^{r) is proportional to the density wj{r), given by the same 
radial form in eq. ([9]) but with a smaller slope parameter b'. This relation underlines the 
close ties between the different potentials. 

To facilitate the calculation of the binding energies, separate bound state solutions were 
computed for V^g^r) and V2g{r) using the fermionic wave functions ipsg^r) ~ Ws{r) and 
'4'2g{r) ~ u'^(r). Further, bound state energies for excited states in V2g{r) were deduced 
by fitting modified Airy functions [9] to the potential. This appears to be a reasonable 
approximation, confirmed by the overall self-consistency achieved. 



The mass of the system is defined by the energy to balance binding 

Mr = -i^sT + ^2," + 2m, , (15) 

where E^'g is the negative binding energy in V^g{r) and Vr!^g{r) (for these potentials only 
the lowest state is discussed here), E2g positive binding energies in V2g(r), and 2mq the 
mass of the elementary fermions (as assumed in the quark model). This shows two types 
of mass generation, binding in the Coulomb like potential V^glr) and dynamical mass 
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generation in V2g{r). The lowest binding energies E^g and i?2^ can be related to the 
average momentum of the bosonic density < q >= J qdq w^{q)/ J dq w'^{q) by the energy- 
momentum relation, which yields -Egg 0* > +-^3g + 2'^g- From this, a further boundary 
condition is obtained 

Mg,s. = <q> +2mg . (16) 

Now, all conditions are specified. However, by calculating relativistic bound states for 
different flavour systems with fermion masses rrig taken from the quark model, a consistent 
solution could be obtained only for light qq systems, as a; (782), for which rriq is assumed 
to be very small. In case of c- or b-quarks, a satisfactory solution could not be found. 

However, an overall internally and externally consistent solution for all systems studied is 
obtained by assuming massless elementary fermions (quantons). The internal consistency 
refers to the mass parameter m, which is in this case the reduced mass due to the binding 
energies in V2g{r) and V3g{r). For the matrix element A^sg (containing two bosons and 
two fermions) the reduced mass m is given by m = whereas for A^2g = 1-^23- 

Therefore, the reduced mass m should be between the limits 

lE2g>m>^\Elg\ or ^ E^g < m < ^ \Elg\ . (17) 

This relation is used as additional constraint. 

The external consistency concerns the description of different flavour systems, which in 
case of massless fermions cannot be distinguished by the fermion masses m^. However, 
in our formalism different flavour systems are described by different slope parameters b 
in eq. (Q, which are related by a vacuum potential sum rule [10]. In this way all flavour 
systems are connected. 

The need for massless fermions has important implications. First, the vacuum of the 
applied theory is the absolute vacuum with average energy E^ac = 0. This is consistent 
with the low energy density of the universe deduced from astrophysical observations. 
Second, in the potential V2g{r) the lowest energy solution is the vacuum and therefore 
Eg = 0. By this, the absolute height of V2g{r) is fixed. Third, by rewriting eq. ( ITTj) in the 
form Ai2g = ■^w'g{q){ip{p')ip{p)}d'^w'g{q)^ one can see that fermion-antifermion pairs can 
be created during the overlap of two fluctuating boson fields. By this mechanism stable 
particles can be created out of the absolute vacuum. 
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In the whole formahsm there are finally five constraints ([7]), f|T4|) . f|T5|) . f|T6|) and f|T7|) . 
by which all five open parameters, two slope (or size) parameters b and b', a shape 
parameter k, the mass parameter m and the coupling constant a, are determined within 
certain ambiguities. In addition, the relation of the slope parameter b of different flavour 
states by a vacuum sum rule [lO] indicates that in principle a complete solution of the 
relativistic bound state problem for all flavour states is achieved. 



An evaluation of the above formalism is discussed for two mesonic systems, u{782) and 
charmonium J/\E'(3097) including excited states. The potentials V^g^r) and V2g{r) have 
been determined by adjusting the open parameters to the flve constraints discussed above. 
Remaining uncertainties have been reduced by adjusting the conflnement potential V2g{r) 
to flt the spectrum of radial excitations, which allows a rather precise determination of 
the mass parameter m. Nevertheless, some ambiguities remain, mainly between k and a. 

Results on the radial dependence of densities and potentials are given in flgs. 1 and 2 by 
using the parameters in table 1. In the upper part the interaction Vy{r) is given by the 
solid line. Compared to the Coulomb potential Vcouii^) = h/r (dot-dashed line) there are 
no divergencies for r — ?• and oo, in agreement with the demand of a flnite theory. 

In the middle part a comparison of the density Wg{r) (dot-dashed line) with the potentials 
V^g{r) (dashed line) and V^g{r) (solid line) is made. We see that condition ((Tj) for the 
vector potential is well fulfllled at larger radii. This indicates that the bosonic wave 
function Ws{r) is well described by the radial form in eq. ([9]) and that also relation ([8]) 
between Ws{r) and w^(r) is correct. 

In the lower part the scalar two-boson density Wg{r) (upper dot-dashed line) is compared 
to the sum of all potentials. Only a considerably smaller slope parameter b' meets the 
boundary condition ( !T^ . showing a good consistency of our description. The resulting 
masses, momenta and mean radius squares are given in table 1 together with the flnal 
parameters. For the determination of a, the angular momentum and spin coupling coeffi- 
cients < qq \L\ 1^ > (for L=0 and 2) are taken into account as well as the corresponding 
components v^{r) of the multipole-expansion of the interaction Vy{r) = J^t^vif^)^ given in 
the present case by f^(r) =< qq \L\ 1^ >^ Vy{r). 

The results for w(782) are given in flg. 1, those for charmonium in flg. 2. Although the 
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Table 1: Results for two mesonic systems, a; (782) and charmonium J /ip including excited 
states, in comparison with the data [3j. Masses are given (in GeV), h and h' in fm, the 
mean radius squares in fm^. For w(782) ^2^=0.51 GeV and E|g=-0.27 GeV, for J/^/^(3097) 
^20=0.34 GeV and Et=-2.76 GeV. 



System 


Ml M| M| M| 


Ml 


^exp ^exp ^exp 


w(782) 

J/^ 
System 


0.78 1.42 1 
3.10 3.69 4 

h V 


.93 

.16 4.58 

K m 


3.2 
12.7 

a 



3 


.782 
.097 


1.42±0.03 

3.686 4.16±0.02 

<rl,> < rl, > 


u;(782) 


0.465 0.220 
0.097 0.050 


1.42 0.25 1.05 
1.39 0.26 0.83 


0.14 
0.10 


0.16 0.04 
0.71 10-2 0.19 10-2 



* aeff = ajdr v^{r)/ J dr Vcouiij) 



binding energies are quite different, in both satisfactory agreement of the various 

quantities is obtained, which fulfil all boundary conditions. It is important to mention that 
in a transformation to momentum space agreement of density and potential is obtained 
only for massless elementary fermions. 

Apart from the requirement of massless fermions, surprising results are obtained also for 
the potential V2g{r), which is shown in fig. 3. This has the same form as the known 
'confinement' potential Vconf = —a/r + 1 ■ r deduced from potential models [IT]. Also the 
comparison with the 'confinement' potential from lattice simulations of Bali et al. [12] 
(solid points with error bars) shows a very similar behaviour. This indicates clearly that 
the confinement of hadrons is a general property of relativistic bound states and can 
be generated in an Abelian or non-Abelian theory. The question of a vector or scalar 
structure of the confinement potential can be studied by looking at the splitting of p- 
wave states in charmonium and bottonium, see ref. [13]. From the existing data neither a 
vector nor a scalar structure is found, supporting the derivative structure of the potential 

As a last point a brief discussion of mass distributions is made (related to the widths and 
decay probabilities), which goes beyond a static description of bound state properties. 
Experimentally, the widths of the mesonic states considered are quite different, 8.4 MeV 
for the first and 90 keV for the second one [3]. The small width of J/ip is explained in 
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the quark model by a different quark; structure involving c-quarks, which cannot decay 
into light quarks. In the present model the width is related to the dynamical structure 
of the system, given by two kinetic energy distributions T2g{q) and Tsg{q). Their Fourier 
transforms T2g{r) and T^gij) are related to the potentials V2g{r) and V^g{r) by the virial 
theorem T2g(r) = V2g{r) fmdT^gir) = i < > {^d'^V^gir) I dr'^+l dV^g{r)/dr). Fromr2g(g) 
and T^g{q) the decay probability of state 1 to state 2 is given by P(f_^2) ~ \^2g{l) ^Igl?)]"^^^ 
and p(f_^2) = \Tlj{q) Tig{q)Y^^ J dQ Yi,{e,(f))Yi,{e,(f)) together with mixed terms. 

Results are given in fig. 4, which show that the overall width of the mass distribution is 
dominated by the two-boson kinetic energy T2g{q) (related to the confinement potential). 
In the Fourier transformation^ this yields a sharp peak with a width of about 0.4 and 0.5 
GeV for the low and high mass system, respectively (solid lines). Differently, T^g{q) yields 
a wide mass distribution, which (including threshold effects in the form used in ref. [14] ) 
is given by dot-dashed lines. 

In both cases the experimental width of the fundamental state is much smaller than 
obtained in our calculations. This can be understood by a rather stable 1~ g.s. with a 
strongly delayed decay. Indeed, a;(782) decays almost exclusively to three pions, the decay 
to two pions is strongly suppressed. 

In the case of J /ip the experimental width of the g.s. is four orders of magnitude smaller 
(~ 90 keV) than the narrow peak in the calculation, indicating a system with a strongly 
delayed decay. The decay probability P(j/^->i^(782)) is less then 10~^° for T2g{q) and ex- 
tremely small for T^g{q). This is mainly due to the large width of T^g{q), which leads 
only to a tiny decay fraction of 4 ■ 10~^ from the region of the 90 keV wide J/ip peak. 
Qualitatively this yields decay widths in agreement with the data, but more detailed work 
is needed to calculate specific decay channels with good precision. Here, only one decay 
should be mentioned, J /ip ^ e^e~ with a decay width of 5.16 keV. By taking the results 
for the electron from ref. [I5] and applying a folding of two electron distributions, a decay 
probability for this channel of 0.8 ■ 10"^ is obtained, corresponding to a width of 4.8 keV, 
which is close to the experimental value. 



^for convergence of the long range potential V2g{r), a Gaussian smoothing is applied. 
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In summary, a solution of the relativistic bound state problem for hadrons has been de- 
rived from a second order Lagrangian, which may be considered as the most natural way 
to overcome the infrared divercency of the SM. The application to qq mesons has led to 
surprising results, a consistent description of different flavour systems by assuming mass- 
less elementary fermions (quantons). This has the consequence that the Higgs- mechanism 
as well as supersymmetric extensions of the SM (in which the infrared divergency prob- 
lem has not been removed) are not needed. Further, the known confinement potential is 
well understood in an Abelian theory (without colour degree of freedom). Our solution is 
self-consistent, without free parameters, and shows a relation to the absolute vacuum. It 
may therefore correspond to a rather fundamental description of hadrons with a simpler 
symmetry structure (less degrees of freedom) than the present (first order) form of the 
SM. 

In a further paper [16j an application of the present formalism to heavy qq bound states 
in the "top"-mass region is discussed with results on the structure of Z{91.2 GeV) and 
the new scalar state at 126 GeV, discovered recently in Atlas and CMS data (tentatively 
interpreted as Higgs-boson). 

The author is very grateful to many colleagues for fruitful discussions, valuable comments 
and the help in formal derivations. Special thanks to P. Decowski and P. Zupranski 
for numerous conversations and encouragements and B. Loiseau for his help with the 
formulation of the Lagrangian and matrix elements. 
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Figure 1: Self-consistent solution for a mesonic system corresponding to uj{19>2). 
Upper part: Interaction Wy{r) in comparison with the Coulomb potential (for a=0.3), 
given by solid and dot-dashed lines, respectively. Middle part: Bosonic density wlij) 
and potential |V3'^(^')| given by the overlapping dot-dashed and solid lines, respectively, 
matched by the condition ([7]); |V3*g(T')| is shown by dashed line. Lower part: Potentials 
V2g{r) (lower dot-dashed line), |V^g('")| (dashed line), |V3'^(r)| (lower solid line) and sum 
(solid line) in comparison with w1{r) (upper dot-dashed line) to meet condition (fUl) . 
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Figure 2: Same as figure 1 for the charmonium system J/'0(3O97). 
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Figure 3: V2g{r) for the two mesonic systems u;(782) and charmonium J/ip, given by 
dashed and dot-dashed hne, respectively, and comparison with the confinement potential 
from lattice gauge calculations [I2] (solid points). 
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Figure 4: Mass distributions for the two systems. The dominant contribution is due to 
the two-boson kinetic energy term T2g{q) (narrow sohd hnes), whereas T^g{q) is given by 
wide dot-dashed distributions (enhanced by a factor 10). The widths are much larger 
than found experimentaUy, indicating a strongly reduced decay probability. 
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